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ABSTRACT 

We present a discussion of possible sources of C, O, and Si in the intergalac- 
tic medium (IGM) using the yields of very massive stars (VMSs) and Type II 
supernovae (SNe II). The chemical evolution of the IGM is considered based 
on analytical phenomenological models of hierarchical structure formation. Two 
regimes are considered: one for gas expulsion by VMSs and SNe II in low-mass 
halos prior to dissociation of H2 molecules and reionization, and the other for 
later SN Il-driven outflows from intermediate-mass halos. We use recent data on 
the abundances of C, O, and Si in the IGM inferred from two UV background 
(UVB) models. We show that the results from a UVB model including quasars 
only cannot be explained by existing stellar models. To account for the results, in 
particular [Si/C]iGM, from a softer UVB model requires VMSs to provide > 15% 
of the C in the IGM. The preferred scenario is that VMSs in low-mass halos 
provided between 15% and 60% of the C in the regime of very high redshift 
[z > 15) and that galactic outflows provided the remainder during later epochs 
(4 < 2; < 6). Thus, there is a large gap in z between metal production by 
VMSs very early in the chemical evolution of the IGM and subsequent contribu- 
tions from galactic outflows. The observational estimate of [0/H]igm implies a 
high [Fe/H]iGM > —3 for all cases considered (including the case of a pure VMS 
source). This raises problems with regard to observations of metal-poor stars in 
the Galaxy. In addition, the [0/Fe] values for SN II models are in conflict with 
stellar observations, which indicate that the calculated average Fe yields of SNe 
II are a factor of ~ 2 too high. These issues and the general problem of relating 
abundances in the IGM to those in metal-poor stars remain to be investigated. 

Subject headings: galaxies: formation — intergalactic medium — nuclear reac- 
tions, nucleosynthesis, abundances 
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1. Introduction 

We present a discussion of possible sources of C, O, and Si in the intergalactic medium 
(IGM) with particular consideration of very massive (> 100 Mq) stars (VMSs) and Type 
II supernovae (SNe II). It is now well established that there is a substantial inventory of 
metals in the general IGM over a wide range of redshift z — 1.5-5.5 (e.g., Songaila 2001; 
Pettini et al. 2003). Detailed studies of the data over z — 1.8-4.1 by Schaye et al. (2003) 
and Aguirre et al. (2004) show that at a given density the abundances of C and Si follow 
a lognormal distribution with substantial scatter. Both the median and the scatter of the 
distribution are dependent on the density but nearly independent of z. Using observations 
at z 2.5 Simcoe, Sargent, & Ranch (2004) also found lognormal distributions for C and 
O. As only the ions C IV, O VI, and Si IV are directly measured, the inferred elemental 
abundances and in particular the relative abundances of C, O, and Si are sensitive to the 
model of UV background (UVB) used (e.g., Giroux & ShuU 1997). We will use the results for 
two common UVB models adopted by the above two groups (see Table 1). For UVB model 
Q including quasars only, Schaye et al. (2003) found the average C abundance in the IGM to 
be [C/H]iGM ~ —2.3 over the density range of the data used (log 5 = —0.5 to 1.8 with S being 
the density relative to the cosmic mean) and up to [C/H]igm ~ —2.0 when extrapolated to 
higher densities. The corresponding [Si/C]iGM is 1-45 (Aguirre et al. 2004). There is no 
evidence for evolution of either [C/H]igm or [Si/C]iGM over z = 1.8-4.1 covered by the data. 
For a similar UVB model, Simcoe et al. (2004) found [C/0]igm = and their lognormal 
distribution corresponds to [C/H]igm = —2.2. For a different UVB model QG including both 
quasars and galaxies, Schaye et al. (2003) and Aguirre et al. (2004) found [C/H]igm = —2.8 
to -2.5 and [Si/C]iGM = 0.74 while Simcoe et al. (2004) found [C/0]igm = -0.5 and a 
median abundance of [C/H] = —3.1, which would correspond to [C/H]igm = —2.5 for a 
lognormal distribution with the same scatter of 0.75 dcx as for model Q. For both models Q 
and QG, the results of the two groups for [C/H]igm are in good agreement. Note that model 
Q gives higher values for [C/H]igm and much higher values for [C/0]igm and [Si/C]iGM than 
model QG (see Table 1). Considering uncertainties in both the measurements and UVB 
models, the average abundances in the IGM arc rather high with [0/H]igm = —2.3 to —2.0. 
The C/0 and Si/C ratios in the IGM do not appear to vary with density (Simcoe et al. 
2004; Aguirre et al. 2004). Both [C/H]igm and [Si/C]iGM appear to be independent of z over 
at least the range oi z — 1.8-4.1. 

An important question is when the metal inventory of the IGM was provided. At present 
we know that this inventory was achieved prior to 2; ~ 4. Metal ejection from both low-mass 
(~ 10^-10^ Mq) halos at early epochs (z > 15) and intermediate-mass (~ 10^-10^° Mq) 
halos at later epochs {z < 15) will be discussed. The stellar sources of these metals are also 
not known but VMSs and SNe II are assumed to be the plausible candidates. We earher 
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proposed a prompt metal inventory from VMSs (e.g., Wasserburg Sz Qian 2000). Using 
estimates of this inventory and yields of Heger & Woosley [2002 (HW)], Oh et al. (2001) 
suggested that VMSs could be the source of reionization (cf., Venkatesan & Truran 2003). 
It will be shown that if the IGM inventory was obtained at early epochs, the required metal 
production rate for VMSs and zero-metallicity SNe II in low-mass halos is much higher than 
the Galactic rate. In contrast, the bulk of the IGM inventory can be provided at later epochs 
by low-metaUicity SNe II in intermediate-mass halos through efficient outflows for a Galactic 
metal production rate. However, in all cases and for all the stellar models considered there is 
a problem in accounting for the elemental ratios inferred for the IGM from UVB model Q. It 
will be shown that the ratios from model QG can be matched with significant contributions 
from VMSs. In particular, the Si/C ratio reported for the IGM appears to require a blend 
of contributions from both VMSs and SNe II if Salpeter initial mass functions (IMFs) are 
assumed for these sources. 

We will not consider the contributions from intermediate-mass (~ 1-8 Mq) stars in the 
following discussion. Such stars are thought to be an important source of C in the Galaxy 
on timescales of > 1 Gyr, thereby contributing to the increase of [C/0] in Galactic stars 
from ~ -0.5 at [0/H] < -1 to ~ at [0/H] ~ (Akerman et al. 2004). Clearly, the 
delayed C contributions from intermediate-mass stars can be ignored for enrichment of the 
IGM at early epochs {z ^ 15 or < 0.26 Gyr since the big bang). The precise epoch where 
intermediate-mass stars contribute to the IGM is not known. However, there is no evidence 
for evolution of either [C/H]igm or [Si/C]iGM over at least the range oi z — 1.8-4.1 (Schaye 
et al. 2003; Aguirre et al. 2004). Therefore, we will only consider VMSs and SNe II as 
the possible sources for metals in the IGM. As intermediate-mass stars produce C but no 
Si, the requirement of VMS contributions to account for the rather high [Si/C]iGM will be 
strengthened by any possible contributions from intermediate-mass stars. 



2. Metal Ejection from Low-Mass Halos 

We first explore the possibility that the metal inventory in the IGM is the result of 
astration (i.e., star formation) in low-mass (~ 10^-10^ M©) halos at 2; ~ 15-30. It is widely 
recognized that the first stars were formed in low-mass halos where gas could cool by H2 
molecules when the virial temperature T^i,. reached ~ 300 K (e.g., Couchman & Rees 1986). 
Haiman, Rees, & Locb (1997) showed that the soft UVB produced by the first massive stars 
could dissociate H2 molecules in low-mass halos universally before the IGM was reionized. 
Dissociation of H2 molecules drastically changes the condition for onset of astration to Tyir ~ 
10^ K as required for cooling by atomic species. This condition then governs astration for 
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all subsequent epochs including that of reionization. Ciardi, Ferrara, Sz Abel (2000) showed 
that universal H2 dissociation has not yet taken place for z ~ 20. While the actual domain 
of astration by H2 cooling remains to be established, we consider the regime of 2; ~ 15-30 
to be within current estimates. Our purpose is to explore if such astration could possibly 
provide the metals in the IGM and what metal production rates are required to achieve this. 

The nature of the first stars made of big bang debris is not known. They may be VMSs 
(see Abel, Bryan, & Norman 2003; Bromm & Larson 2004 for reviews). On the other hand, 
the discovery of a low-mass star with [Fe/H] = —5.3 (Christlieb et al. 2002) suggests that 
low-mass stars may have formed very early (e.g., Schneider ct al. 2003) and presumably a 
broader spectrum of normal stars including SN 11 progenitors might be expected for the first 
stars. In either case, the occurrence of a VMS or SN II in a low-mass halo would disrupt the 
baryonic gas in the halo and disperse it into the IGM on some scale. To estimate the range 
of low-mass halos with gas expulsion, we use models of hierarchical structure formation in 
the standard cold dark matter cosmology as discussed in Barkana & Loeb (2001). 

A critical parameter for halo evolution is the virial temperature T^ij. of the gas. For 
z > 1, Tvir is related to the halo mass M (mostly in dark matter) as 



where is the mean atomic weight and n = 1.22 or 0.6 for a neutral or ionized gas, respec- 
tively, with a primordial composition of H and He. We assume that astration starts when 
the gas reaches a threshold virial temperature Tvir,o for cooling. For H2 cooling we will take 
Tvir.o = 300 K {11 = 1.22). Another critical parameter is the gas binding energy £'bi,gas- For 



We assume that all the gas would be expelled from the halo by an SN II or VMS with 
an explosion energy E^^p if £'bi,gas < -E'exp (e.g., Bromm, Yoshida, & Hernquist 2003) and 
retained for -E'bi,gas > E^^p. For illustration, we take E^^p — 10^^ erg for SNe II and 4 x 10^^ 
erg for VMSs (HW). 

For a halo associated with an na density fluctuation (an na halo), the growth of its mass 
M{z) as a function of z is prescribed by the Press-Schechter formahsm (Press & Schechter 
1974; see Fig. 6 in Barkana & Loeb 2001). Astration starts in a 3a halo when the halo 

mass reaches Mq = 3.6 x 10^ M© at z = 26.9. The gas in the halo would be expelled by 
an SN 11 until the halo mass grows to Mu = 7.0 x 10^ Mq at z = 18.6 or by a VMS until 
Mbi = 7.1 X 10^ Mq at ^ = 15.3. In general, at a given z an nga halo would reach M{z) = Mq 




2/3 



(1) 



z > 1, 
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to start astration and an nuc halo would reach M{z) — Mbi to greatly hinder gas expulsion. 
At any given z all na halos with uq < n < corresponding to Mq < M < have 
started astration but will be disrupted by an SN 11 or VMS. The relevant range of halos at 
z = 20 corresponds to 2.3 < n < 3.2 (5.5 x W Mq < M < 6.7 x 10^ Mq) for SNe 11 or 
2.3 < n < 3.8 (5.5 x W Mq < M < 6.1 x W Mq) for VMSs. For a given z, the fraction 
F{Mq < M < M\yi\z) of all matter in na halos with no < n < uu is 

F(Mo <M <Mu\z) = ^J " exp (^'y) (3) 

The function F(Mo < M < Mu\z) for T^ir,o = 300 K and Eu,gas = 10^^ erg is very close 
to that for the same rvir,o but -Ebi.gas = 4 x 10^^ erg as halos with -E'bi,gas > 10^^ erg are 
extremely rare for 2; > 15. A fraction ~ 0.05% to 5% of all matter is in low-mass halos with 
gas expulsion at 2; = 30 to 15. 

We now estimate enrichment of the IGM by gas expulsion from low-mass halos using O 
as a representative element. We consider a large reference region of the universe and treat 
it as a closed homogeneous system. We assume that the expelled gas is immediately mixed 
with the entire IGM of the reference region. In this case, the number ratio (0/H)igm of O 
to H atoms in the IGM as a function of time t (since the big bang) is determined by 

^J2I3iSU ^ Po,„/(H)„,„ (4) 

n=no 

where Po,no- is the O production rate in an na halo and the sum extends over all the na halos 
with no < n < in the reference region. We explicitly assume that Po^na is proportional 
to the number (H)^^. of H atoms in the gas of the na halo, 

Po,na = Ao(0/H)q(H)„,. (5) 

For z ^ 1, the majority of the H atoms in the reference region reside in the IGM and the 
majority of the H atoms in a halo reside in the gas. This gives 

^ Ao(0/H)qF(Mo <M< Mbil^). (6) 

The constant parameter Aq is related to the frequency of occurrence {R) for the O source 
in a halo of mass M as 



«.AoXg^^l44Ao^^J^3^). (7) 



6 



where Xq — 9.6 x 10~^ is the solar mass fraction of O, fb — 0.15 is the baryonic mass fraction 
for the halo, and (Yq) is the average O yield of the source. Equation (6) can be integrated 
to give 



where t(z) = 0.538[10/(1 + 2;)]'^/^ Gyr for the adopted cosmology and t' = t{z'). The above 
treatment is readily extended to mixtures from different sources. 

The result for [0/H]igm = logZ^^ as a function of z is shown in Figure 1 for Aq = 
0.1 Gyr~^ using Tvir,o = 300 K for Mq and -Ebi.gas = 10^^ erg for Mbi (solid curve). Replacing 
the mass range Mq < M < Mbi with M > 10^ Mq gives essentially the same result (dot- 
dashed curve). The rate Aq = 0.1 Gyr~^ is representative of O production by SNe II in 
the Galaxy [the solar O abundance is produced by SNe II over ~ 10 Gyr]. For this rate 
[O/HjiGM = —3.5 at 2; = 15. So SNe II in low-mass halos cannot provide the IGM inventory 
with a nominal abundance of [0/H]igm — —2.3 unless the relevant production rate is 16 
times higher than the Galactic rate. The requirement of Aq = 1.6 Gyr~^ would also apply 
to any other sources (including VMSs) if they were to produce the IGM inventory at early 
epochs. Simulations by Norman, O'Shea, and Paschos (2004) appear to suggest that VMSs 
could not enrich the IGM to the level of [0/H]igm = —2.3. However, this result was based on 
putting one VMS in each halo with M > 5 x 10^ Mq that was found in a simulation volume 
of 1 Mpc^ at z = 15. The number of such halos found in the simulation volume may not 
be statistically representative (M. L. Norman, personal communication). Future numerical 
studies on early enrichment of the IGM by VMSs or any other sources should use larger 
simulation volumes to check the required metal production rate against the value given here. 



We next consider contributions to the IGM from just galactic outflows at 2; < 15. We 
assume that only normal stars including SNe II progenitors are formed and astration starts 
when Tvir.o = 10^ K {/i — 0.6) is reached for cooling by atomic species. The halo mass 
at the onset of astration ranges from Mq = 4.6 x 10''' Mq for 2; = 15 to 2.6 x 10^ M© for 
2; = 4. For these masses -Bbi.gas > 10^^ erg so an SN II cannot expel all the gas from the halo. 
Subsequent to the onset of astration, we assume that a fraction e of the metals produced by 
SNe II is lost from the halo and injected broadly into the IGM until the halo reaches a cutoff 
mass Ml. While the dependence of e on the halo mass may be complicated, we will treat 
it as a constant to determine how efficient mass loss must be to provide the IGM inventory 





3. 



Galactic Outflows 
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prior to z — A. Applying the same formalism leading to equation (8), we have 

rt{z) 

Zjp^ = eAo/ F{Mo < M < Mi\z')dt', (9) 

where ti5 = 0.26 Gyr is the time corresponding to 2; = 15. Using Aq = 0.1 Gyr~^ and 
e = 0.1, we show [0/H]igm = logZ^^^ as a function of z for Mi = 10^ and 10^° M© (short 
and long dashed curves), respectively, in Figure 1. These Mi values are consistent with 
observations of outflows from Lyman break and dwarf galaxies (e.g., Pettini et al. 2001; 
Martin, Kobulnicky, & Heckman 2002). Values of Mi > 10^° Mq give almost the same result 
as Ml = 10^° M©. Figure 1 gives [0/H]igm = -2.9 at z = 4 for Mi = 10^° Mq. Thus to 
achieve [0/H]igm = — 2.3at2; = 4by just galactic outflows requires a rather high but still 
reasonable e = 0.4 for a Galactic metal production rate of Aq = 0.1 Gyr~^ or an equivalent 
eAo = 0.04 Gyr~^. Obviously, smaller values of eAo are required if galactic outflows are to 
provide some signiflcant portions but not all of the IGM inventory (see §4) . 



4. Elemental Ratios in the IGM and Implications for the Sources 

The results for in §§2 and 3 can be extended to other elements using the relevant 

production rates. For two elements E and E', 



Ae (1e) /X| 



Ae' {Ye') 



(10) 



where 

(Ye) = .Z J I ■ (11) 

j^^ (p{m)am 

In the above equation m is the mass in units of Mq for the source with mi and m„ being the 
lower and upper limits, iE(m) is the yield of element E as a function of m, and 4>{m) = is 
the IMF. We assume a Salpeter IMF (/? = 2.35) for SNe II. As there is no a priori knowledge 
on the mass distribution of VMSs, a wide range of IMFs with /3 = 2.35, 10, and 20 are 
explored. This allows us to examine a wide range of abundance ratios of Si and Fe relative 
to C and O resulting from VMS sources. The average yields (1e) for C, O, Si, and Fe and 
the corresponding [C/0], [Si/C], and [0/Fe] are given in Table 2 for a number of models 
of VMSs [Umeda & Nomoto 2002 (UN); HW] and SNe 11 [Woosley & Weaver 1995 (WW); 
UN; Chieffi & Limongi 2004 (CL)] for the j3 values considered. Clearly, the elemental ratio 
E/E' for a source only depends on those stars with finite yields of E and E'. For this reason, 
we give the limiting masses for metal production as mi and m„ in Table 2. Although some 
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VMSs or SNe II may form outside this mass range, they do not affect the calculation of the 
elemental ratios. 

Table 2 shows that VMSs have a narrow range of [C/0] — —0.71 to —0.50 for a wide 
range of IMFs but their [Si/C] varies from 1.18-1.30 for a Salpeter IMF to 0.61-0.84 for an 
IMF with P — 20 that very strongly favors the lower masses. This is because the C and O 
yields are nearly constant for all contributing VMSs while the Si yield increases significantly 
(the Fe yield increases steeply) with mass (UN; HW). In contrast, models of zero-metallicity 
SNc 11 (WW; UN; CL) give narrow ranges of [C/0] = -0.23 to 0.02 and [Si/C] = 0.13-0.17. 
For SNe II with very low to solar metalhcities, WW give [C/0] = -0.42 and [Si/C] = 0.42- 
0.50 while CL give [C/0] = -0.2 and [Si/C] = 0.18-0.28. We consider that SNe II in the 
regime of low metalhcities are relevant to the problem at hand. In all cases, [Si/C] is high 
for VMSs but low for SNe II with zero or low metalhcities. It is the high value of [Si/CjioM 
that indicates significant VMS contributions as noted by Aguirre et al. (2004) (see also 
Schaerer 2002). As VMSs have [Si/C] > [Si/C]iGM while SNe II have [Si/C] < [Si/C]iGM, 
this supports the existence of VMS contributions but was not discussed in the recent critique 
of VMS models by Tumlinson, Venkatesan, and ShuU (2004). 

We first compare the elemental ratios from stellar models with [C/0]igm = and 
[Si/C]iGM = 1-45 inferred from UVB model Q. Neither VMSs nor SNe II of any mctallicity 
can account for both [C/0]igm and [Si/C]iGM- As [C/0] is close to [C/0]igm for SNe II 
of zero metallicity (or all mctallicitics according to CL) and [Si/C] is close to [Si/C]iGM for 
VMSs with a Salpeter IMF, one may consider a mixture from these two sources. However, 
to approximately match [C/0]igm would require that almost all of the C come from SNe II 
and this would give a [Si/C] for the mixture essentially equal to the SN II value, which is off 
by > 1 dex from the IGM value. Similarly, matching [Si/C]iGM would require that almost all 
of the C come from VMSs, thus giving a very low [C/O] for the mixture. So in all cases and 
for all the stellar models there is a problem in accounting for the elemental ratios inferred 
from UVB model Q. 

We now consider UVB model QG with [C/0]igm = -0.5 and [Si/C]iGM = 0.74. These 
ratios cannot be matched by VMSs with a Salpeter IMF, for which [Si/C] > 1.18. A match 

can be obtained for VMSs if an IMF very strongly favoring the lower masses is used (see 
models 2B, IC, and 2C in Table 2). It is also possible to obtain a match with both [C/0]igm 
and [Si/C]iGM for UVB model QG if we consider mixtures of VMS and SN II contributions. 
As the largest [Si/C] = 1.3 for VMSs (model lA) and the largest [Si/C] = 0.42 for SNe 
II of low metallicity (model 6), to match [Si/C]iGM = 0.74 requires a fraction of at least 
jVMS _ g "i^g ^jjg Q come from VMSs. Below we discuss four representative cases, all of 
which match [C/0]igm and [Si/C]iGM for UVB model QG to within 0.1 dex (see Table 3). 
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In case I the IGM inventory was produced by just VMSs (IMF with /3 = 20) in low- 
mass halos at early epochs {z > 15). A high metal production rate of Aq''^^ = 1.6 Gyr"-*^ is 
required to achieve [0/H]igm = —2.3 (§2). This corresponds to a frequency of occurrence 
Ryus = 4(M/105 Mq) Gyr-i (eq. [7] for (Y^^^) = 54.3 Mq) so that every halo of ~ 10^ M© 
would have formed ~ 1 VMS by 2; = 15 (tis = 0.26 Gyr). This rate is usually considered 
the upper limit for VMS formation as the UV irradiation from a VMS suppresses further 
astration in its hosting halo (e.g., Couchman & Rees 1986). Case I does not allow any 
significant contributions to the IGM from galactic outflows at later epochs {z < 15). This 
appears to be in conflict with observations of outflows from Lyman break and dwarf galaxies 
(e.g., Pettini et al. 2001; Martin et al. 2002). It is also at odds with the attribution of the 
bulk of the IGM inventory to galactic outflows as advocated by many previous studies (e.g., 
Madau, Ferrara, & Rees 2001; Scannapieco, Ferrara, & Madau 2002; Thackcr, Scannapieco, 
& Davis 2002). In addition, Bromm & Loeb (2003) suggested that VMSs may form only 
for metallicities below [C/H] = —3.5 ±0.1 and [0/H] = —3.0 ± 0.2. However, for case 1 
VMSs are required to provide [C/H] = —2.9 and [0/H] = —2.3, which are much beyond the 
suggested termination point for their formation. 

In case II the sources for metals in the IGM are a combination of VMSs and zero- 
metallicity SNe II in low-mass halos at early epochs. A single Salpeter IMF (/3 = 2.35) 
is used for the mass range covering both VMSs and SNe II. To achieve [0/H]igm = —2.3 
requires 

^VMS ^ 0.9 Gyr-^ and A^^" = 0.7 Gyr^^ The latter is 7 times greater than the 
Galactic value. Thus case II requires high metal production rates for both VMSs and zero- 
metaUicity SNe II. As in case I, case II does not allow signiflcant contributions from galactic 
outflows and requires VMSs to provide enrichments beyond the suggested termination point 
for their formation. 

For cases III and IV we consider a combination of early contributions from VMSs (IMF 
with P — 10 for case HI and Salpeter IMF for case IV) and later galactic outflows driven by 
low-metallicity SNc II. Case HI is chosen to match [Si/C]iGM = 0.74 exactly. The required 
rate of Aq^^^ = 1.0 Gyr~^ is high but the required galactic outflows with an efficiency of 
e = 0.1 for a Galactic rate of Ag^" = 0.1 Gyr'^ (for Mi = 10^^ Mq, see §3) arc quite 
reasonable. Case HI also requires VMSs to provide [C/H] = —3.0 and [0/H] = —2.5, 
which are in excess of the termination metallicities for their formation as suggested by 
Bromm & Loeb (2003). Case IV was calculated using Jq^^ — 0.15, which requires VMSs to 
provide metals only up to the suggested termination point, and hence, a much smaller rate 
of Aq^^ — 0.3 Gyr~^. The required efficiency of galactic outflows is e = 0.3 for a Galactic 
Aq^^^. If we instead had required an exact match to [Si/CjicM for case IV, then /q^^ — 0.23 
with no major difference in the results. 
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The results in Table 3 show that VMSs must be significant contributors to the IGM 
based on the [Si/C]iGM value and that either a pure VMS source with an extreme IMF 
or a blend of VMS and SN II sources could provide the observed metal inventory in the 
IGM. However, case IV appears to be the best scenario if we require (1) major contributions 
from galactic outflows, (2) a Galactic Aq^'^^ (3) a nonextreme IMF for VMSs, and (4) strict 
termination metallicities for VMS formation as suggested by Bromm & Loeb (2003). In 
case IV VMSs in low-mass halos provided 15% of the C, 20% of the O, and 51% of the Si 
prior to H2 dissociation and reionization, and the rest of the IGM inventory was provided 
by later efficient outflows driven by SNe II in intermediate-mass halos. This suggests that 
VMSs must be substantial contributors to the IGM with high frequencies of occurrences at 
very high z. Venkatesan & Truran (2003) calculated that VMSs would provide 10 H-ionizing 
photons per baryon in order to reach [0/H]igm = —2.5. So the VMSs in case IV would 
provide ^ 3 H-ionizing photons per baryon in reaching [0/H] = —3.0 (using C gives the 
same result. Note that the value of 0.35 H-ionizing photons per baryon given in Venkatesan 
& Truran 2003 assumes a termination metallicity of [0/H] ~ —4.0 for VMS formation). This 
may be consistent with the "weak" version of the VMS scenario discussed in Tumhnson et 
al. (2004). 

So far our discussion has focused on C, O, and Si as these elements are measured in 
the IGM. The abundances of other elements such as Fe can also be calculated for the four 
cases discussed above. The O and Fe yields in Table 2 give [0/Fe] = 0.71 and 0.52 for 
cases I and III, respectively, and [0/Fe] = 0.12 for cases II and IV (see Table 3). For 
[0/H]iGM = —2.3, this imphes [Fe/H]iGM = —3.0 and —2.8 for cases I and III, respectively, 
and [Fe/H]iGM = —2.4 for cases II and IV. Observations of metal-poor stars in the Galaxy 
(e.g., Cayrel et al. 2004) show that [0/Fe] > 0.5 for [Fe/H] ~ -3.0 to -2.4. Cases I and 
III are consistent with these observations. However, the [0/Fe] value for cases II and IV 
is too low compared with those observed in metal-poor stars. This may be a potentially 
important problem for these two cases, especially for case IV, which satisfies all the other 
criteria discussed above. The low [0/Fe] value for cases II and IV is mostly due to the 
low values of [0/Fe] ~ 0.2 given by SN II models. We note that the Fe yields of SNe 
II were not calculated a priori but generated by employing largely artificial mass cuts in 
all models. Stellar observations (e.g., Nisscn et al. 2002) show that [0/Fe] ~ 0.5 for the 
regime of —2.4 < [Fc/H] < —1.5, where SNe II are clearly the dominant source for both O 
and Fe. This indicate that the calculated average Fe yields of SNe II are a factor of ~ 2 
too high. Much work is required to resolve the problem of Fe yields and the associated 
uncertainties in [0/Fe] for SN II models. We also note that our next best case HI gives 
an adequate [0/Fe] but requires VMSs to provide [C/H] and [0/H] somewhat beyond the 
termination metallicities suggested by Bromm & Loeb (2003). If the criteria for termination 
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of VMS formation could be relaxed, then case III may be a viable alternative to case IV in 
accounting for the IGM data and stellar observations. 

Clearly, the results presented here are subject to uncertainties in theoretical models of 
stellar yields and in IGM measurements and analyses. Possible uncertainties are unlikely 
to explain the large discrepancies between stellar models and the elemental ratios from 

UVB model Q but may affect some of our results regarding model QG. As discussed above, 
uncertainties in Fe yields of SN II models affect the [0/Fe] values for cases II and IV. In 
addition, the required rates quoted above (see Table 3) correspond to a nominal abundance 
of [0/H]iGM = —2.3 and would change if different values were used. For this nominal value 
all four cases considered (including case I for a pure VMS source) imply [Fe/H]iGM > —3.0 
(using [0/Fe] = 0.5 instead of 0.12 for cases II and IV does not change the result). This is 
high considering a number of ultra-metal-poor stars in the Galaxy have —4 < [Fe/H] < —3 
(e.g., Christlieb 2003). One possible explanation would be that those ultra- metal-poor stars 
formed in small halos in under-enriched regions and were later accreted into the Galaxy. 
This is consistent with the large scatter for the lognormal distribution of metal abundances 
in the IGM and with the standard model of hierarchical structure formation. The above 
issues and the general problem of relating abundances in the IGM to those in metal-poor 
stars remain to be addressed. 

In summary, we find that VMSs must be significant contributors to the IGM. Our 
preferred scenario is that VMSs provided between 15% (case IV) and 60% (case III) of 
the C (and corresponding levels of other metals) in the IGM by disrupting the gas in low- 
mass (~ 10^-10^ Mq) halos at very high redshift {z > 15). Outflows driven by SNe II in 
intermediate-mass (~ 10^-10^° Mq) halos provided the remainder (the bulk in case IV) of 
the metals during later epochs (4 < 2; < 6, see Fig. 1). Thus, there is a large gap in z 
between metal production by VMSs very early in the chemical evolution of the IGM and 
subsequent contributions from galactic outflows. The metals from these two different sources 
appear to be well mixed locally but have a wide range in net abundances as established by 
Schaye et al. (2003), Aguirre et al. (2004), and Simcoe et al. (2004). 

We thank the referee, Aparna Venkatesan, for a thorough report and many helpful 
suggestions that greatly improve the paper. This work was supported in part by DOE grants 
DE-FG02-87ER40328 (Y. Z. Q.) and DE-FG03-88ER13851 (G. J. W.), Caltech Division 
Contribution 9010(1115). 
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Table 1. IGM Data^ 



UVB Model 


[0/H]iGM 


[C/0]iGM 


[Si/C]iGM 


Q 


-2.3 to -2.0 





1.45 


QG 


-2.3 to -2.0 


-0.5 


0.74 



1C/0]iGM from Simcoe et al. 2004, [Si/C]iGM from 
Aguirre et al. 2004, and [0/H]igm from Simcoe et 
al. 2004 or converted from [C/H]igm of Schaye et 
al. 2003 using [C/0]iG.\i- Solar abundances adopted: 
log £0(0) = log(C/H)0 + 12 = 8.52, log €0(0) = 8.83 
(Grevesse & Sauval 1998), and loge0(Si) — 7.55 (An- 
ders & Grevesse 1989). 



Tabic 2. VMS and SN II Yields and Yield Ratios 





Metallicity 








{Y^) {Mq) 










Model 


{Z/Zq) 


(mi,m„) 




C 





Si 


Fe 


[C/O]^ 


[Si/C]^ 


[0/Fe]*' 


VMSs 






















lA (UN) 





(150, 270) 


2.35 


3.67 


51.0 


18.4 


6.85 


-0.71 


1.30 


0.00 


2A (HW) 





(140, 260) 


2.35 


4.34 


44.2 


16.6 


6.21 


-0.57 


1.18 


-0.02 


IB (UN) 





(150, 270) 


10 


4.56 


52.4 


12.7 


2.99 


-0.62 


1.05 


0.38 


2B (HW) 





(140, 260) 


10 


4.99 


46.4 


9.24 


0.67 


-0.53 


0.87 


0.97 


IC (UN) 





(150, 270) 


20 


5.41 


54.3 


9.41 


1.42 


-0.57 


0.84 


0.71 


2C (HW) 





(140, 260) 


20 


5.51 


47.2 


5.59 


0.046 


-0.50 


0.61 


2.14 


SNe n 






















3 (WW) 





(12, 30) 


2.35 


0.12 


0.30 


0.039 


0.099 


0.02 


0.13 


-0.38 


4 (UN) 





(13, 30) 


2.35 


0.28 


1.29 


0.098 


0.083 


-0.23 


0.15 


0.32 


5 (CL) 





(13, 35) 


2.35 


0.31 


1.19 


0.11 


0.10 


-0.15 


0.17 


0.19 


6 (WW) 


0.01 


(12, 40) 


2.35 


0.20 


1.43 


0.13 


0.13 


-0.42 


0.42 


0.17 


7(CL) 


0.005 


(13, 35) 


2.35 


0.34 


1.49 


0.13 


0.10 


-0.21 


0.18 


0.28 


8 (WW) 


1 


(11, 40) 


2.35 


0.19 


1.37 


0.15 


0.12 


-0.42 


0.50 


0.20 


9 (CL) 


1 


(15, 35) 


2.35 


0.36 


1.54 


0.17 


0.14 


-0.20 


0.28 


0.19 



"^Calculated using the same solar abundances as in Table 1. 

''Calculated using X§ = 9.6 x 10"^ and X® = 1.3 x 10"^ (Anders & Grevesse 1989). 
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Table 3. Representative Cases to Match the IGM Data from UVB Model QG 



Case 


VMS 
Model 


SN II 
Model 


fVMSa 

/c 


[C/0] 


[Si/C] 


fVMSa 

Jo 


[0/Fe] 


^VMS ^SNII 

(Gyr-i)b (Gyr-i)b 




I 


IC 




1 


-0.57 


0.84 


1 


0.71 


1.6 






II 


lA 


4 


0.28 


-0.42 


0.82 


0.54 


0.12 


0.9 


0.7 




III 


2B 


6 


0.60 


-0.49 


0.74 


0.66 


0.52 


1.0 


0.1 


0.1 


IV 


2A 


6 


0.15 


-0.45 


0.65 


0.20 


0.12 


0.3 


0.1 


0.3 



^Specified by the same Salpeter IMF governing both VMSs and SNe II for case II, chosen 
to match [Si/C] = —0.74 exactly for case III, and chosen to obtain [0/H] = —3.0 from VMSs 
for case IV. 

''Rates and efficiencies required to provide [0/H] — —2.3. The required rate Aq^^ for 
VMSs in all cases and the required rate A^'^" for zero-metaUicity SNe II in case II are given 
assuming complete disruption of the gas by these sources in low-mass halos. For galactic 
outflows in cases III and IV, the required outflow efficiency e is given assuming a Galactic 
Aq^" (an equivalent eAg^^^ also satisfies the requirement). 
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Fig. 1. — Evolution of [0/H]igm as a function of z for gas expulsion from low- mass halos with 
astration controlled by H2 cooling (solid curve for contributing halos with Mq < M < Mbi 
determined by T^i^i^ = 300 K for onset of astration and -E'bi,gas = 10^^ erg for onset of 
gas retention, dot-dashed curve for contributing halos with M > 10^ Mq) or outflows from 
intermediate-mass halos subsequent to H2 dissociation (Tvir,o = 10"^ K, short and long-dashed 
curves for cutoff masses of 10^ and 10^° Mq, respectively). The vertical line at 2; = 15 divides 
the two regimes. For both regimes the reference production rate Aq = 0.1 Gyr~^ corresponds 
to the Galactic rate for SNe II. 
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